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In this paper we address the problem of consistently construct Langevin equations to describe 
fluctuations in non-linear systems. Detailed balance severely restricts the choice of the random 
force, but we prove that this property together with the macroscopic knowledge of the system is 
not enough to determine all the properties of the random force. If the cause of the fluctuations 
is weakly coupled to the fluctuating variable, then the statistical properties of the random force 
can be completely specified. For variables odd under time-reversal, microscopic reversibility and 
weak coupling impose symmetry relations on the variable-dependent Onsager coefficients. We then 
' analyze the fiuctuations in two cases: Brownian motion in position space and an asymmetric diode, 

, for which the analysis based in the master equation approach is known. We find that, to the order 

of validity of the Langevin equation proposed here, the phenomenological theory is in agreement 
with the results predicted by more microscopic models. 

o 



in 



I. INTRODUCTION 

[ A widely used method to study fluctuations is the so-called Langevin approach, introduced by P. Langevin ||l| as a 
way to study Brownian motion |^,||. In his treatment, the interaction force between the Brownian particle and the 
] ^ bath particles is split into two contributions. The first one corresponds to the damping due to the frictional force 
. exerted by the bath. The second one arises from the thermal motion of the bath molecules. The Brownian particle 
Qv^ ■ then experiences a large number of collisions per unit of time with the bath molecules, which give rise to a rapidly 
: varying contribution, responsible for the observed erratic motion of the particle. Then, the equation for the velocity 
d " of a Brownian particle of unit mass can be written as 



T' u{t) = -^uit) + Fit) (1.1) 

-a ■ 
a ■ 
o 

Q The first term of the right hand side is the friction force. This force is a linear function of the variable u, ^ being the 
i;^ friction coefficient. The second term stands for the Langevin complementary force or random force, whose nature is 
unknown a priori. Beyond Langevin's original treatment and with the aim of determining u{t), several assumptions 
have been made on the nature of the random force. The most basic one is that F{t), being an unpredictable rapidly 
\ varying function of time, is taken as a stochastic process, often assumed to be Gaussian and white with zero mean 
1^. The second moment of F{t) is determined by means of the fluctuation-dissipation theorem, which relates the 
strength of the random force to the dissipation, i.e. the friction coefficient, ^. This crucial point ensures that the 
law of equipartition of energy is satisfied by the system under the effect of F{t). With all these assumptions, the 
statistical properties of the random variable u{t) can be obtained 

Let us consider here a generic macroscopic variable, X{t), whose relaxation is phenomenologically described by the 
equation 

X{t) = -A{X{t)) (1.2) 



>< 



where A{X{t)) is the flux of the variable which can be a non-linear function of X{t). In this paper, we will be 
interested in the dynamics of the fluctuations of that variable. The fluctuating counterpart of X{t) will be denoted 
by x{t), the former being obtained from the latter by some averaging procedure. It has been then proposed that x{t) 
satisfies the equation 

±{t) = -A{x{t)) + F{t) (1.3) 
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obtained by replacing X{t) by x{t) in the phenomenological equation (1.2) and adding a Langevin force, assumed to 
be Gaussian and white. This approach has been used for systems in which A(x) is a linear function of the variable x 
[Q-^. However, when A{x) is non-linear, the addition of a random force F{t) to eq. (l.S), as a procedure to describe 
the random process x{t), has often been questioned. We want to emphasize two main criticisms [Q: 

1. The Gaussian nature of the random process F{t) is postulated but a more general character of the process 
could in principle be equally plausible. 



2. Averaging eq. (1^) with respect to the realizations of the F{t) given an initial condition xq^ identifying X{t) 



with {x{t))x„ and comparing with the phenomenological equation (1.2), one finds that 



{A{^m.„ + A{X{t)) (1.4) 

if A(x) is not a linear function of x. It is the n ar gued that this difference arises from the fact that the use of 
A{x) also in the equation for the fluctuations ( |l.3| ) is not justified because A(x) is phenomenological, that is, it 
is known except a function of the order of the strength of the fluctuations|^. 

In view of these criticisms, it is often concluded that the Langevin approach is not valid for the case of internal 
(thermal) noise in non-linear systems and that a more microscopic description is necessary to derive the proper 
equation for the fluctuations 

Recently, however, some attention has been payed to the properties of the Langevin random force in the case of 
thermal noise for a physical system bearing a non- linear A{x). Let us write A{x) under the form 

A{x) = ~i5{x)^\^P,{x) (1.5) 
where P^jx) is the equilibrium probability density of the variable x, supposed known, and /3(a;) is a given function of 



Eq. (1.5) emphasizes the fact that the macroscopic flux is a function of the thermodynamic force 



T{x) ^ ^\xiP,{x) (1.6) 

such that the flux vanishes when the thermodynamic force also vanishes. Near equilibrium, tends to a constant and 
A(x) takes the form of a linear function of the thermodynamic force T ^ in the spirit of Onsager's theory of irreversible 
processes. Thus, iOx) is the corresponding variable-dependent Onsager coefficient. P. Mazur and D. Bedeaux have 
proved that the underlying microscopic reversibility when applied to the Langevin equation, severely restricts the 
properties of the Langevin complementary force, contrary to what is stated in the first remark. Indeed, if x(t) in eq. 
(|1.3D is a variable even or odd under time-reversal and if F{t) is independent of the variable x{t), it is found that 
[^: 1) F{t) must be a Gaussian process and 2) P{x) must be a constant. Note that if Pe{x) is not Gaussian, then 
A(x) is non-linear even if /3 is a constant and however, F{t) still must be a Gaussian process. Furthermore, another 
conclusion that can be drawn from this result is that if f3{x) is not a constant, F{t) cannot be independent of x(t). 
This case will be referred to as genuinely non-linear. 

A simple way to construct a random force depending on the variable x(t) is to assume that ^0 11 

F{t) = £{x{t))L{t) (1.7) 



often referred to as multiplicative noise. Here, L{t) is a stationary white process, independent of x{t) and c{x) is a 
given non-vanishing function of x. 



^See also the comments in ref. M about this point. 
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In general, when proposing a Langevin equation to study thermal fluctuations in a given system, one would wish 
that the properties of the random force would follow from the phenomenological knowledge of the system. It is clear, 
as claimed in ref. that if a detailed microscopic description is available, there are ways to derive the equation 
for the fluctuations of a given variable without ambiguity. However, we are interested in the situations where this 
information is lacking, although it will be assumed that there exis ts a n underlying microscopic description in terms 
of Hamilton or Schrodinger equations, and that a relation like eq. (^]^) is known. The aim of this paper is to use the 
information at hand to consistently determine a Langevin equation for x{t) in the genuinely non-linear case. To this 
end and in the same spirit as in refs. |^- we will use generally valid properties such as microscopic reversibility 
to fix the nature of the random force. Furthermore, we will discuss the range of applicability of the theory developed. 
To that purpose, we will consider that the equation for the evolution of x{t) is given again by the phenomenological 
law ( pT^ ) plus two additional contributions. The first one is a variable-dependent random force, F{x{t),t), defined by 

F{x{t),t)= I dx5{x-x{t))L{x,t) (1.8) 



where L{x,t) is a stationary and white random process, depending on a parameter x but independent of random 
process x{t), which is characterized by its cumulants |^ 

< L{x,t) >c = (1.9) 

< L{xi,ti)L{x2,t2) >c =C2{xi,X2)S{ti-t2) (1.10) 
< L{xi,ti)L{x2,t2) . ■ ■ L{Xn, tn) >c = Cn{xi,X2, ■ ■ ■ , Xn)S{tn - ti)S(tn - ^2) • • ■ S{tn ~ tn-l) 

(1.11) 



Note that while L{x,t) is independent of the process x{t), the statistical properties of F{x{t),t) will depend on the 
evolution equation for x{t). 

The second contribution is a function T(x), standing for a possible modification of the phenomenological flux A{x) 
due to fluctuations. The proposed Langevin equation then reads 

x{t) ^ -A{x(t)) + t{x(t)) + F{x{t),t) (1.12) 



A priori, both r(a:) and the set of functions Ci{{xi}) are undetermined. We will specialize in the simplest, although 
very important, case in which the random process L{x,t) depends only on one single function c{x), all the cumulats 



present in eqs. (1.9 1.11) being related with c{x). This case will be later identified with the widely used multipHcative 
noise 0,0. The general case will be analyzed elsewhere. The presence of the unknown function r(a;) in eq. (|1.12| ) 
as well as the treatment to be developed are the main differences with respect to the analysis given in refs. [pj-p^ 1[ , 
permiting a larger degree of freedom with important consequences. 



The paper is organized as follows. In the next section we will obtain the master equation for the one-dimensional 



process x(t) described by eq. (1.12) with multiplicative noise. In section 3 we apply detailed balance to the master 
equation to show that L{x,t) needs to be Gaussian, which constitutes an alternative derivation of the results of ref. 
[ pT[ . In addition, we will also obtain a differential equation involving c(x) and f (a::), as well as the symmetries imposed 
to these functions, reminiscent of the Onsager symmetry relations. In our context, such an equation serves to fix 
neither c{x) nor r{x) but it is a mere relationship between them. This sugges ts th at detailed balance alone is not 
enough to univocally determine the properties of the additional terms in eq. ( |l.l2| ). In section 4, we will analyze 
Brownian motion in position space to illustrate the main ideas exposed so far. This particular example will serve us 
to ilustrate that if the mechanisms causing the fluctuations are weakly coupled to the fluctuating variable, then we 
can univocally determine c{x) and r(a;). This will be called the weak coupling assumption. In the same section, we 
also apply these ideas to Alkemade's diode |ljjl^, often used to show the failure of the phenomenological treatment 
of the fluctuations. This last example will permit us to analyze the phenomenological theory developed along the 
paper to the light of a more microscopic description and, therefore, determine its range of applicability. Section 5 is 
devoted to the conclusions and to a brief summary of the ideas sketched in this paper. 
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II. THE MASTER EQUATION 



Let us assume that the dynamics of the one-dimensional variable x{t), even or odd under time-reversal, is described 



phenomenologically, that is, in the absence of fluctuations, by eq. (1.2). Note that this statement constitutes a 
definition of the phenomenological equation. For the fluctuations we have the Langevin equation given in eq. ( 1.12| ) 
which we rewrite as 

xit) = -B{x{t))+F{{x{t),t) (2.1) 



with 

B{x) = ~(3ix)^\nP,ix)-t{x) (2.2) 

where the form ( pT^ ) for A{x) is used. Note that the " renormalized" flux B{x) is not completely determined since it 
depends on the unknown r(a;). We will define an acceptable B{x) as such that, after fixing r(a;), still retains some 
dependence on the macroscopically relevant Onsager coefficient /3(x). 



The random force F{x(t),t), as it stands in eq. (1^), cannot be univocally interpreted unless a rule to compute 
averages of the form < g{x{t))L{x,t) >, g{x) being an arbitrary function, is provided In the literature, this 
ambiguity is often referred to as the Ito-Stratonovich dilemma. For causality reasons, we demand [plfol that the state 



of the system at t is not correlated with the random force at the same time, that is 

< F{x(t), t) dx< 6{x - x{t)) >< L{x, t) >= (2.3) 



which in fact prescribes Ito's rule. However, from the mathematical point of view, the Ito-Stratonovich dilemma is 
immaterial in our case since T{x), and the properties of L{x,t) are a priori unknowns that will be determined later 
on, on the basis of physical arguments and in agreement with the chosen interpretation rule. 

For simplicity's sake, in the calculation which follows we will assume Stratonovich's rule. Then, for clarity, we 
introduce r(x), F{x(t),t), L{x,t) and the cumulants {c„(a;i, . . . ,x„)}, as the quantities consistent with Stratonovich's 
rule while r(a;), F{x{t),t), L{x,t) and {cn{xi, . . . ,Xn)} will be those consistent with Ito's rule. Since both describe 
the same physical situation, we must have that the equations 

xit) - -A{x{t)) + T{x) + F(a;(t),t)|stratonovich (2-4) 

x{t) = -A{x{t)) + t{x) + F{x{t), t) (2.5) 

ltd 

are equivalent. 

Assuming thus Stratonovich's rule, we introduce the multiplicative noise as the particular case in which the set of 



cumulants in eqs. (1.9-1.11) takes the form 

Cn{xi,X2, ■ ■ . ,Xn) = a„c(a;i)c(x2) . . . c(.T„) (2.6) 

Q!„ being constants. This particular choice of the random process L{x, t) makes it equivalent to c{x)L{t) as in eq. 



(1.7), with L{t) a s tati onary white random process whose cumulants are the set of constants {a^}. Note that the 



choice done in eq. (2.6) does not imply a Gaussian nature of the random process L{x,t) which, in turn, is going to 



be Gaussian if = for i > 3. 

To obtain the master equation for the probability, we will essentially follow ref . |^ . Let us first consider the density 
distribution 
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p[x, t) = 5{x — x{t)) 



(2.7) 



where x{t) is the solution of eq.( p.4[ ) for a particular realization of L{x, t) and a given initial condition Xq = x{t = 0) . 
Then, the conditional probability, P{x^t\xQ), for the random variable x{t) to be at the point x at t, given that it was 
at the point xq a,t t = Q, is related to p{x, t) by 

P{x,t\xo) =< p{x,t)>^^ (2.8) 



where, again, the averages are over all the realizations of L{x, t) and the subindex Xq indicates that the same initial 
condition xq is considered. Clearly, 



F(x, t — 0\xo) — 6{x — Xq) 



(2.9) 



In appendix A we show that P{x, t\xo) satisfies 



d 

—P{x,t\xo) 



o oo 1 / r) \ " 

^BW + i:(-l).'-a„ -cW 



n=2 



P{x,t\xa) 



(2.10) 



which is the Kramers-Moyal expansion of the master equation 
d 



-g^P{x, t\xo) = / dx' [w{x' x)P{x', t\xo) ~ w{x x')P{x, t\xo)] = / dx' W{x\x')P{x' ,t\xo) (2.11) 



Here, w{x x') stands for the transition probability from the state x to x', and the operator >V(x|x') is defined as 

>V(2;|x') = w{x' — > a;) — 6{x — x) I dx" w{x x") (2-12) 



Eq. ( 2.10 ) can be expressed in a more compact form by introducing the operator C{x): 



£(.).-B(.) + E(-l)";^-(a^^(-) 

n— 2 ^ 



(2.13) 



Then 



d_ 
dt 



Pix,t\xo) = C{x)P{x,t\xo) 



(2.14) 



We thus identify 



W{x\x') ^ £{x) 6{x - x') 



(2.15) 



Furthermore, due to the Markovian nature of the process under discussion, we can write 

P{x,t)^ / dx' P{x,t\x',t')P{x',t') 



(2.16) 



where P{x, t ) is t he one-point probability density. Setting x' = xq and t' = in this last equation, multiplying both 
sides of eq. ( 2.14 ) by P(xo,0) and integrating over xq one gets the evolution equation for P{x,t) 

^^P{x,t)^C{x)P(x,t) (2.17) 
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III. DETAILED BALANCE 



A consequence of the time-reversibility of the microscopic equations of motion in a closed system in equilibrium is 
t he so -called detailed balance p4|. When the dynamics of the variable x{t) is described by a master equation like eq. 



(2.11), detailed balance is the condition imposing that for each pair of states x and x' the transitions must balance: 

w{x' x)Pe{x') — w{ex ex')Pe{ex) (3-1) 



with e = 1 for even and e ~ — 1 for odd variables under time-reversal. Furthermore, the equilibrium probability 
distribution satisfies Pe{x) = Pe{ex). 



Detailed balance as sho wn in eq. (3J) leads to a similar relationship for the operator yV(a;|a;') defined in the 
previous section (eq. ( 2.15| )) |^ 



W{x\x')Pe{x') = Wiex'\ex)Peix) (3.2) 

which, in turn, leads to the relationship 

C{x) Pe{x) ^{x) = Peix)C^{ex)i;{x) (3.3) 

where i'{x) is an arbitrary function. The operator jf^ix) is the adjoint of jC{x) in the usual sense and is given by 



(3.4) 



Even variables 



Le t us first dis cuss the case of even variables and thus replace e by 1 in eq. (3.3). Making use now of eqs. ( 2.13| ) 
and (3.4) in eq. (3.3) we get 



= P.(.)|-i.(.)|- + |:l«„(c(.)A)"|^(.; 



(3.5) 



We now multiply both sides of this last equation by c{x) obtaining 



I VB{x) + f^i^ir^anV^cix) I Pe{x)i;{x) = Pe{x) I -B{x)V + f2 ^anc{x)V" I V(x) 

L n=2 ' J L n=2 ' J 



(3.6) 



where the we have introduced the operator 2? as 



V = c{x) 



dx 



(3.7) 



6 



Expressing the left hand side of eq. (3.6) in terms of the independent operators V^, one has 

n=2 ' m=0 ^ ^ 



(3.8) 



Since il}{x) is arbitrary, we can now equate the coefficients of the terms proportional to 'D"^'ip{x) in both sides of this 
last equation. From the coefficient of Tf'ijjlx) = we obtain 



oo 

[DB{x)Pe{x)] + 5](-l)"-a„ [V^c{x)Pe{x)] = 



(3.9) 



Dividing this equation by c{x) and using eq. ( |3.7| ), one gets the stationarity condition C{x)Pc{x) = as follows from 
eq. ( ^.17 ). Similarly, from the coeffcient of 'Dil>{x) we get 



B{x)Pe{x) + ^(-l)"-a„n [2?"-ic(a:)Pe(x)] = -B{x)Pe{x) 

n=2 



(3.10) 



Finally, the coefficient of T)'^4'{x) for m > 2 gives 

V ^(-l)"a„ f " ) [P"-™c(x)Pe(x)] - —^a^c{x)P,{x)- for m > 2 
■■^ n! \ m J <- ■'to! 



(3.11) 



This last expression stands for an homogeneous system of equations for the set of constants at for i > 2 whose 
coefficients are functions of x. Note that ay, can only appear in the equation for to = 2 in view of the lower limit of 
the sum on the left hand side of eq. ( 3.11 ). However, for to = 2 the term containing a2 in the sum is exactly canceled 
by the term on the right hand side, so that 02 is not determined by this system of equations. In a ppendix B we proof 
that the only acceptable solution of eqs. ( 3.11 ) is 0;^ = for i > 3. Therefore, L{x,t) in eq. (2^) must be Gaussian. 
This important result was already found in ref. mU by using a different procedure. 



Thus, eqs. ( ^.g| ) and ( |3.10 ) now take respectively the form 



|-{5(.)P.(.) + ia. 

2B{x)Pe{x)+a2 



cix)—c{x)Pe{x) 

d 

c{x)—c{x)Pe{x) 



= 





(3.12) 
(3.13) 



according to the definition of V. Notice that eq. ( p.l2| ) can be obtained from eq. (3.13) by differentiating with respect 
to X. While the former is the stationarity condition which states that the divergence of the probability flux must be 
zero in equilibrium, the latter says that the probability flux itself must also vanish. 



Odd variables 



The treatment of the case of odd variables cannot be carried out following the same procedure as done for even 
variables for a general non- vanishing c{x). It is, however, of particular importance the situation in which c{x) verifies 
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c{x) = c{—x) since we will prove that L{x,t) is Gaussian if and only if c{x) is an even function of its argumentg. 
Furthermore, we will also determine the corresponding symmetry properties of B(x) and discuss the consequences 
that arise from these facts. 



Let us replace e by —1 in eq. (3.3) and assume that c{x) is even. We then get 



= Peix) |B(-a:)^ + f;(-irl«„ (^c(x)^y'l^(x) (3.14) 



Using again the operator V and equating the coefficients of the functions 'D'"'^ip{x) at both sides of eq. (3.14), we 
arrive at 

oc 

[DB{x)P,{x)] + V(-l)"^a„ [D'^c{x)P,{x)] = (3.15) 

n=2 

for the coefRcient of 'D'^ijj{x). From the term proportional to 'Dip{x) we get 

oo ^ 

B{x)Pe{x) + V(-l)"-a„ n [V"-^c{x)Pe{x)] = B{~x)Pe{x) (3.16) 

n=2 

The coefficient of T>™iIj{x) for m > 2 gives 

f2 ^(-l)"«n f " ) [D--"'c{x)Pe{x)] = (-l)™i^a„c(x)Pe(x); for m > 2 (3.17) 



In appendix B we also proof that the only acceptable solution of this last equation is — for i > 3, so that L{x^ t) 
is again a Gaussian process. 

The reverse implication is not difficult to proof. Let us assume that L{x, t) is a Gaussian process or, what is the 



) is i 

In 



same, that cti = for i > 3. From cq. ( p.3| ) with e = —1, making use of eqs. ( p. 13 ) and ( |3.4D we obtain 



= Pe{x) I B{-x)— + -a2 ( c(-x)— ) ) V(x) (3.18) 



where now c{x) is a general non- vanishing fimction. Equating the coefhcients of the terms proportional to d^ijjix) / dx"^ 
we get 

]^a2C^{x)Pe{x) = ]^a2C^{-x)P,{x) (3.19) 



^This was implicitly assumed in eq. (2.7) of ref. [ pi] . 
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and since c{x) is assumed non-vanishing, then c{x) = c{—x). Therefore, if detailed balance is satisfied and x(t) is an 
odd variable under time-reversal, L{x,t) is Gaussian if and only if c{x) is an even function. 



Making now use of the Gaussian nature of L{x,t), eqs. ( 3.15 ) and ( 3.16| ) become 



§-[B[x)P.(x) + \a, 
{B{x) ^ B{-x))Pe{x) + a2 



d 

c{x)—c{x)Pe{x) 
c{x)—c{x)Pe(x) 



= 

= 



(3.20) 
(3.21) 



Again, eq. (3.20) has the structure of the divergence of a flux equal to zero. Since this flux itself is zero in equilibrium, 
we get 



B{x) = ~B{-x) 



(3.22) 



Thus, under these conditions eqs. (3.2C) and (3.21) reduce to eqs. (3.12) and (3.13) and, therefore, under these 
circumstances the subsequent discussion will be valid for even as well as for odd variables. 



While eq. ( |3^ ) follows fr om detailed balance, we will prove in the next section that if weak coupling is satisfied, 
this previous symmetry implies in turn 



T{x) = ~V{~x) 
Pix) = p{-x) 



(3.23) 
(3.24) 



The Fokker-Planck equation 



Let us consider again eq. ( |3.13 ). In view of eq. (2^) it can be written as 



( 



dx 



2dx 2 



(3.25) 



Note that the factor 02/2 may be absorbed i nto c ^jx). In what follows, we will, without lost of generality, choose the 
scale of L{x,t) such that a2 = 2. Then, eq. (3.25) turns into 



^(^) _ ^(^)) ^ inP,(a;) = r(x) - ~c'{x) 



(3.26) 



With the aim of discussing the original Langevin eq. ( 1.12 ) or, equivalently eq. (2_^), we have to relate the previous 
results with those corresponding to ito's prescription. It is easy to show although we will not do it here, that for 
r(a;), c{x) and L{x,t), one can find that i) L{x,t) is also Gaussian with c.2{x,x') = c{x)c{x'), ii) T{x) and c.{x) are 
related by the equation 



(3.27) 



and, for odd variables, iii) the same symmetries as in eqs. ( |3. 22 - 3. 24 ) hold 



Eq. ( 3.2g ), bearing two unknown functions r(a;) and c{x), indicates that for our Langevin equation, microscopic 
reversibility alone cannot univocally determine the properties of the random force F{x(t),t). Rather, it relates the 
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functions T{x) and c{x) with the thermodynamic force and the phenomenological coefficient (3{x). In the next section 
we will introduce a plausible way to determine these two functions based on physical grounds. 



Eq. ( 3.27 ), the fact that the random process L[x,t) must be Gaussian as well as the symmetries of f (a;), c{x) and 
j3{x) for odd variables, constitute the main results of this section. Some consequences follow. 

Firstly, in view of the fact that L{x, t) must be a Gaussian process, then eq. ( |1.12| ) is equivalent to a Fokker-Planck 
equation of the form 

^P{x, t)^^ (b{x) P{x, t) + ^^{x) P{x, t)] (3.28) 



Secondly, integrating eq. (3.27) we get 

c^x) 



n 



Pe{x) 



■f3{x) 



Peix) 



rfx'Pe(x'){^/3(x')-f(x') 



f3{x) - A(x) 



(3.29) 



where is an integration constant. In this last expression we have gathered in A[x) all the terms that explicitly 
depend on the equilibrium probability distribution Pe{x). Note that, in general, the properties of the random force 
are then dependent on the Onsager coefRcient (3{x) as well as on the equilibrium probability distribution. 



Inserting this last expression for c^{x) into eq. ( |3.28 ), after some algebra, one gets 



dt 



P{x,t) 



dx 



iPix) - A(x)) 



-P{x,t) — \nP,{x) + 



dx 



P{x,t) 



(3.30) 



Note that this equation depends on the still unknown function f (x) in view of eq. ( [3.29| ). Nevertheless, independently 
of that function, the equilibrium probability distribution Pe{x) is always a solution of this equation. 

Finally, after remark 2 and due to the fact that r(a;) is not necessarily zero for a linear A{x), we have no objective 
reason at this point to consider the theory for linear A(x) as substantially different from that of the genuinely non-linear 
case. 



IV. TWO EXAMPLES 



In this section we will analyze two non-trivial examples for which a more detailed analysis exists. The first one 
is the diffusion of a Brownian particle with a position-dependent friction coefRcient in a general potential field. The 
second one is a diode with a non-symmetric non-linear characteristic function / = I(V) 



Diffusion of Brownian particles 



Here, we will study the case in which the variable x{t) is the position of a Brownian particle in a one-dimensional 
system in equilibrium with a reservoir that keeps constant the temperature T. The particle is sensitive to a potential 
field V{x). Thus, the equilibrium distribution function is given by 

Fe(x) = -ie-^(-)/'=^ (4.1) 



where J\f is the normalization constant and k is Boltzmann's constant. The phenomenological equation in this case is 
the relationship between the velocity of the particle and the force acting on it 
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xit) 



1 



V{X) ^ ~AiX) 



(4.2) 



£,{X) being the friction coefScient which is assumed to be position-dependent. Physical examples could be either 
motion in an inhomogeneous medium, or an analogy with the three-dimensional situation where hydrodynamic inter- 
actions with other particles or with the walls of the container causes the friction coefficient of a Brownian particle to 
change from one point to another. From these two equations we can rewrite A{X) as 



A{X) 



i{x) dX 



^{X) dX 



(4.3) 



which permits us to identify (3( X) 
equation of the form of eq. (^.l[) 



kT/£^{X). To describe the fluctuations, let us replace X by x and rewrite an 



i(t) = 



kT 



^{x{t)) dx 



\nPe{x)+V{x{t)) 



+ F{x{t),t) 



(4.4) 



where the term between brackets is B{x{t)). If detailed balance has to be satisfied, the corresponding Fokker-Planck 
equation then reads 



d ^, , d f kT 



(4.5) 



Note that this equation does not coincide with the well known Smoluchowski equation that describes the diffusion of 
Brownian particles with position-dependent friction coefRcient (see, for instance, refs. [p|JT^) unless A(a;) identically 
vanishes. 

At this point we will further assume that the properties of the random force need to he independent of Pe (x) (This 
statement will be referred to as "weak coupling" from now on). If this is the case, from eq. (3.27) we get that, in 
general. 



i?{x) 
f(x) 



P{x) 
ax 



(4.6) 
(4.7) 



so that the random force is completely characterized. To get a more intuitive picture of the physical nature of this 
statement, in the case of Brownian motion, we may say that weak coupling is equivalent to state that the random 
force should be a property of the bath system and then it must be independent of the potential force, —dV{x)/dx, 
externally applied to the Brownian particle. Then, if eq. (4.7) is satisfied, A(a;) identically vanishes and the usual 
diffusion equation is recovered. 



Note that in the case of odd variables, since we have assumed that c(x) is even, eqs. ( |4.6D and (4.7) leads to the 
symmetry properties shown in eqs. (3.23) and (3.24). 



We can then conclude that our Langevin equation (1.12) is completely determined when detailed balance and weak 
coupling are imposed, leading us to the Gaussian nature of L{x,t) and to eqs. (4.6) and (4.7). We have also found 
that such a Langevin equation is equivalent to the Fokker-Planck equation 



d_ 

dt 



P{x,t) 



d 



P{x,t) 



d_ 

dx 



InPe(x) 



d_ 

dx 



P{x,t) 



(4.8) 



for the probability distribution. 
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Diode 



We will next consider the problem proposed in refs. [iJJT^, where an extensive analysis on more microscopic grounds 
can be found. One is now interested in the study of charge fluctuations in a condenser of capacity C in parallel with a 
vacuum diode which is in thermal equilibrium with a reservoir at a temperature T. The diode is constructed by facing 
two electrodes of metals with different work function^ for the electrons, Wi and W2 with Wi > W2. The system is 
described in the scheme of fig.l. It is important to note that we have explicitly included in the flgure the contact 
potential barrier, equal to the difference in electrochemical potential per unit charge, when one electron passes from 
one metal to the other. This contact potential is AW/e, where AW = Wi — W2 and e is electron's elementary charge. 
For the system of the figure we have the I{V) characteristics 

I(V) = lo (e"^/'^^ - 1) (4.9) 



The relation between the intensity and the voltage difference v between the diode plates is obtained by using in eq. 



(4.9) the fact that V — AW/e + w as it follows from the figure. It then reads 



One way to proceed is to connect a condenser of capacity C in parallel with the system of the figure 1 to study 
charge fiuctuations in the condenser originated in the diode. In that case, the equilibrium distribution for the charge 
in the condenser follows by standard equilibrium statistical mechanics. Effectively, the energy of the condenser when 
storing a charge q is purely electrostatic and equal to E = g^/2C. Then, the equilibrium distribution is given by 

Pe{q) - e'^ (4.11) 



With the equilibrium distribution (4.11) and the phenomenological law (4.9), one can obtain the equation for the 
fluctuations following the same procedure as in the previous example. The system studied in refs. fl^ , [l^ is, however, 
slightly different and corresponds to that shown in fig. 2. Note that since the system has no contacts, the plates of the 
condenser need to be of different metals. This crucial point leads us to the equilibrium distribution function for the 
charge in the condenser. One considers again the energy of the condenser when a charge q is stored. To discharge 
the condenser, on one hand, the electrons loose electrostatic energy in traveling from metal 2 (where the voltage is 
according to fig. 2) to metal 1 (where the voltage is v = q/C). On the other hand, there is an additional lost of 
energy due to the fact that the chemical potential of the electrons in metal 2 is higher (lower work function) than in 
metal 1 (higher work function). The energy stored in the condenser is thus E = (j^/2C + AWq/e and the equilibrium 
distribution reads 

Pe{q) ^ e ^I'Tc (4.12) 



The equilibrium charge go of the condenser is thus 



AWC 

90 = (4.13) 



These results obtained using equilibrium statistical mechanics are also found in ref. p^ l from the master equation 
proposed to describe the system under discussion. Moreover, the fact that the condenser in equilibrium bears a charge 



The work functions are defined as the work needed to extract an electron from the metal 
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go could be interpreted as a violation of the second law of thermodynamics. It is clear from the fact that the two 
plates of the condenser are of different metals that, in equilibrium, no work can be obtained by connecting them with 
a wire since the equilibrium voltage vq = qo/C = —AW/e would exactly be compensated by a contact voltage barrier 
due to the difference in electrochemical potentials. 



Consequently, in order to study charge fluctuations in this system, we have to use eq.( 1.12| ) together with eq.( l.lO ). 
In view of fig. 2, note that the charge in the cond enser decreases when / is positive. Thus, the phcnomenological 
equation for this system is obtained from eq. (4.10) and reads 



q = -I{v) - /o (l - e^(9-9«)/^-^'^) 



(4.14) 



where eq. ( 4.13 ) has b een used. Finally, with the aim of comparing with the results obtained in ref. ||13[, we rewrite 
eqs. ( 1.12| ) and (4.14) in terms of the dimensionless variable x = —e^^'^{q — qo)/e, where e = e^/kTC is a small 
parameter that is related to the amplitude of the fluctuations^. We get 



(4.15) 



and 



x = ] 



(e-^'''^ - l) 



(4.16) 



where j = /q e^^^/e. As in the previous example, from these two equations we obtain the phcnomenological coefficient 



m = ^ (i - e-^"^) 



(4.17) 



Finally, making use of eqs. (4.15) and (4.17) in eq. (4.8) we arrive at 



xP{x,t) + —P{x,t) 
ox 



(4.18) 



To end this section, we compare eq. ( 4.18 ) with the expansion of the master equation for this model. We thus 
expand our Fokker-Planck equation in powers of e^/^. We obtain 



a' 



dx^ 



l_iei/2^+±e:j2 p(^^^)_^ 0(^3/2) 



(4.19) 



On the other hand, eq.(45) of ref. [n3| has the form 



''This parameter is related to the inverse of the size of the system, in this case the capacity of the condenser C [Q| . Physically, 
e measures the ratio between the difference in electrostatic energy between two consecutive electron jumps and kT. The jump 
probability of a second electron is strongly affected by the jump of the first one if e ~ 1. 
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^P{x, t) = ^(x- \e''\x^ - 1) + \e{x^ - 3x) ) t) 



dx^ 

g3 



l-^e'/^x+^eix'^l)]Pix,t) 



ylex)pix,t) + ^^e^,Pix,t) + Oie^^^) 



(4.20) 



Comparing these last two equations we see that they agree up to order e^/^ while up to first order, the equation for the 
fluctuations is not even a Fokker-Planck equation but higher order derivatives have appeared. We then conclude firstly 
that to the order of validity of the Fokker-Planck equation itself, (e^/^), the fluctuations are correctly described by the 
phenomenological theory developed along this paper without ambiguity. Secondly, up to this order, the coefficients 
of the Fokker-Planck equation are functions of the variable x, corresponding to a so-called non-linear Fokker-Planck 
equation. The possibility of correctly set non-linear Fokker-Plank equations from a phenomenological theory is the 
main result of this paper. As we have seen, however, the validity of the Fokker-Planck equation to describe fluctuations 
in non-linear systems is restricted to small fluctuations. 



V. CONCLUSIONS 



In this paper we have tried to answer the question about whether it is possible or not to write a Langevin equation 
to describe the equilibrium fluctuations of a given macroscopic variab le X with only a phenomenological knowledge 



of the system. To this end, we have proposed the Langevin equation (1.12) containing two terms accounting for the 
existence of fluctuations. In the first place, we have a function T(x) which plays the role of a modification of the 
phenomenological law at the scale of the fluctuations. This answers one of the major criticisms to the use of Langevin 
equation to describe fluctuations in non- linear systems |1. In the second place, we have introduced a causal random 



force as given in eq. (1.8). We have then analyzed the particular case in which the cumulants of L{x,t) depend on 
a single function c{x), which is equivalent to a multiplicative noise, which is the simplest case of variable-dependent 
random force. Then, the function c{x) can be interpreted as an amplitude of the random force that can vary from 
point to point. We have proved that detailed balance leads, firstly, to the result t hat L(x,t) must be a Gaussian 



process. One important consequence of this result is that our Langevin equation ( 1.12 ) is equivalent to a Fokker- 
Planck equation. Secondly, we have obtained the relationship between the functions T{x) and c{x) given in eq. (3.27). 
While these are general results for even variables, we have seen that they are only satisfied for odd variables under 
certain symmetry conditions. In fact, in this last situation we have proved that L{x,t) is a Gaussian process if and 
only if c(x) is an even function, which also implies that B{x) must be odd. The symmetry of this "renormalized" 
flux leads in turn to the fact that r{x) has to be odd and the Onsager coefficient /3{x) has to be an even function. 
This is an important condition on the phenomenological coefficient that arises from microscopic reversibility and weak 
coupling, which has the same origin as the symmetry of the Onsager coefficients for the crossed terms in the linear 
case 

Moreover, we have seen that the phenomenological knowledge of the system together with detailed balance are 
not enough to determine the properties of the random force. The existence of this ambiguity indicates that, at this 
level, different equations for the fluctuations could be proposed, all satisfying detailed balance and having the same 
equilibrium probability distribution, but leading to different dynamics for the fluctuations [ p^ . We have then shown 
that, if the internal mechanism that causes the random force is weakly coupled with the variable that fluctuates, as it 
has clearly been shown in the case of Brownian motion, detailed balance together with this weak coupling assumption 
suffices to completely determine the Langevin equation. 

In our analysis we have imposed weak coupling only at a later stage of the discussion in order to completely specify 
the properties of the random force, after having applied detailed balance. It is important to note, however, that 
if we had imposed detailed balance together with weak coupling from the beginning, we would have seen that the 
Gaussian nature of L{x,t) follows both for variables even and odd under time-reversal, that the "renormalized" flux 
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is authomatically acceptable, and that the symmetry property of /3{x) is a natural consequence of both principle^. 



The results derived from the previous analysis have bee n ap plied to two non-trivial examples. To the range of 
validity of the form assumed for the random force in eq. (1^), the results coincide with more detailed models for 
the dynamics of the fluctuations in those systems. Taking into account that for a general process the validity of 
the Fokker-Planck equation as an approximation of the master equation is restricted to small fluctuations , in the 
second example of section 4 we have seen that, to the range of validity of the Fokker-Planck equation itself, our 
approach is correct and furnishes a description of the fluctuations in a genuinely non-linear system. Therefore, all the 
results found here by means of a phenomenological analysis are consistent with those obtained from the expansion of 
the master equation. 



Let us analyze the case /3 constant in more detail. From eqs. 



and (|4.7D, we get that <? = (3 (also constant) and 



r = 0, in agreement with ref. |^ which, for linear Langevin equations, reduces to the standard results. It is important 
to realize that detailed balance alone e nsure s, in view of eq. ( |3.27 ), that < T{x) >— 0, but still r(a;) and c{x) can be 
non-constant functions related by eq. ( ^.27 ). Therefore, the standard theory for Langevin equations, either linear or 
with constant /?, relies on an implicit weak coupling assumption^. 

On the ba sis o f causality and stationarity, the second moment of the random force in a non-linear Langevin equation 
of the type (^.l|), has been proved to satisfy 



{F{xit),t)F{x{t'),t')) = {L{x{t),t)L{xit'),t'))^2{xit)B{xmS{t~t') 



(5.1) 



Inserting the definition of B{x) in the right hand side of this equation, after partial integration we get 



{x{t)B{xm = {fiixm + (x{t) —pix) 



{x{t)f{xm 



(5.2) 



x—x{t) I 



If detailed balance applies, one can make use of eq. 
cancel and we finally obtain ||l^ 



(4.7). Then, the second and third terms in this last equation 



{F{x{t),t)F{x{t'), t')) = 2{p{x{t))) 5{t ~ t') 



(5.3) 



Again, if /3 is constant, only now we recover the fluctuation-dissipat ion theorem as a c onse quence of both detailed 
balance and weak coupling. In the genuinely non-linear case, eq. (5.3) as well as eq.(4.6) are reminiscent of the 
fluctuation-dissipation theorem existing for linear systems. In particular, eq. (4.6) relates the amplitude of the 
random force, c{x), to the phenomenological coefficient I3{x) at every point. 



Finally, one can calculate the average of the flux in equilibrium from eq. 
(|3.27D it follows that 



(1.12). Effectively, from causality and eq. 



f3{x{t)) — InPe(x) 
ax 



x—x{t) 



+ 

ax 



(5.4) 



x—x(t) I 



as expected. Notice that the underlined term is the modiflcation of the phenomenological law due to the fluctuations 
and ensures that < i: > is zero in equilibrium, so that there is no violation of the second law as it is the case in the 
so-called " Brillouin paradox" Q . 



^This can be proved by treating in eq. (3.3) both ipi^) s^nd Pe{x) as arbitrary functions and following the same reasoning as 
in appendix B but for a general c(x). 

^This implicit weak coupling assumption is formulated by demanding that the random force F(x{t),t) is independent of the 
variable x{t). 
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APPENDIX A: DERIVATION OF THE MASTER EQUATION 



The density distribution given in eq. (2.7) satisfies a continuity equation 



d , . d . , . , , 
— p(a;,t) = -—x{x,t) p{x,t) 



(A.l) 



In view of eq. ( ^.7| ), and due to the properties of the (5-function, to compute x{x, t) we can use eq. (2.4) replacing the 
random variable x{t) by the field variable x, giving 



^p{x, t)^~^ i-B{x) + L{x, t)) p{x, t) 
ot ox 



(A.2) 



where use has been made of eq. (1.8). Eq. (A.2) can be formally solved, obtaining 



p{x, t) = { exp I dt' ^ iB{x) — L{x, t)) ^ S{x — xq) 



(A.3) 



where use has been made of p{x, t — 0) — 6{x — xq). Note that the term between curly brackets is independent of the 
initial condition a;o since x is not a random variable but a field point and the process L{x,t) has been assumed to be 
independent of x{t), s o th at it is also independent of x{t = 0) = xq. Thus, we can average both members of this last 



equation and use eq. (2.8) to also derive the formal solution for the conditional probability density 

P{x, t\xo)^ ^exp dt' ^ {B{x) - L{x, t))^ P{x, t = 0|a;o) (A.4) 



The average is performed over all the realizations of L{x, t) and it is independent of xq. This last expression can also 
be written in terms of the cumulants 

P{x, t\xo) = exp 1 ^ ^ dt' ^ [B{x) ~ L{x, t)]^ ^ | P(x, t = 0\xo) (A.5) 



For n — I one has 







* d 
dt' — [Bix) 



L{x,t)] 







* d 
dt' ^ [B{x)- 



< L{x,t) >, 



(A.6) 



For n> 2 one has 




L{x,t)]] ) ^(-l)«<a„(^c(:E)) (A.7) 
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whe re in deriving the last Hne use has been made of eqs.(l.ll) and ( p.6[ ). Therefore, substituting these results into 
eq. ( |A.5| ) one gets 



P{x, t\xo) — cxp < t 



Bix) + J2{~ir-aJ-cix) 

n=2 ' ^ 



dx 



Pix,t = 0\xo) 



(A., 



Differentiating with respect to time, one arrives at eq. (2.10) 



APPENDIX B: PROOF OF a/ = FOR J > 3. 



In this appendix we will show the possible solutions of the system of equations (3.11) and (3.17) for the parameters 
ai, i > 2. Once removed a2 (recall that the coefficient of the only term involving this constant is zero), the system is 
a homogeneous set of equations for for i > 3 with x-dependent coefficients. Its solution is = for i > 3 if the 
determinant is non-zero and, therefore, L(x, t) is a Gaussian process. However, if its determinant is zero, then other 
solutions with some ai ^ Q may exist. In o rder t o find the expression of the determinant, let us separately write the 
equations for m even and odd. From ( [3.1l| ) and ( 3.17 ) it follows that 



{2k) 



n=2k+2 



= 2k 



1 + e 



{2k- 



'—y{x)a2k+i+ E ^ ( 2fc + 1 ) = rn = 2k + l 

n=2fe+2 ■ ^ 



(B.l) 
(B.2) 



with k = 1,2, .... We have defined the function y{x) = c{x)Pe{x) for the ease of notation. 



Let us first consider the case of even variables. We thus replace e by 1 in eq. (B.2). Making now use of this 
last equa tion, we can express the coefficient of a2k+i in terms of all higher ones. If one then proceeds to substitute 
eqs.(B.2) systematically for all k in eqs. (B.l), one arrives at a set of equations only for the even coefficients. Moreover, 
this new set of equations is triangular and it is then sufficient to calculate the terms of the diagonal in order to know 
its determinant. One can see that the determinant is zero only if 



y{x) {-{Vy{x)f + y{x)Vy{x)) - 



(B.3) 



which implies that some ai may be non-zero if y{x) satisfies the non-linear differential equation 



■c\x) 



c{x)y{x) 



dc dy 

dx dx 



c^{x)y{x) 



fy 

dx'^ 



= 



(B.4) 



whose solutions are 



y{x) = c{x)Pe{x) = V + ji I Pe{x)dx 



(B.5) 



V and /i being integration constants. This last equation determines c{x), since Pe{x) has been considered as given. 
Note that in this particular case, c{x) is independent of the phenomenological coefficient I3{x). Furthermore, these 
solutions are eigenfunctions of the differential operator T), 



Vy{x) = fj.y{x) 



(B.l 
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In order to determine the set of cumulants tti, we will use this last property in eqs. (B.9) and (3.10). We get 

VBix)Peix) 



2B{x)Pe{x) 



n^2 



n=2 



(B.7) 
(B.8) 



Let us introduce the function 



n^2 



(B.9) 



terms of which equations (B/7) and (B^) become 

VB{x)Pe{x) 

y{x) 

2B{x)P,{x) 



(B.IO) 
(B.ll) 



Thus, eliminating B{x) between eqs. ( B.IO ) and ( B.ll ), we obtain a differential equation for i/'(^) 

d/i /i 



(B.12) 



which indicates that ■0 is a quadratic function of the parameter /i. This implies that all the cumulants of L(x, t) should 
be zero ex cept for t he second one. Therefore, we have proved that also in this case L{x,t) is Gaussian. Moreover, 
from eqs. ( B.9| -B.ll) we arrive at the corresponding Langevin equation 



dx c{x) 
It ~ ~Y' 



a2[i + L{x,t) 



(B.13) 



Now, in view of eq. (B.5), we see that eq. ( B.13| ) is independent of the phcnomenological coefficient P{x), corresponding 
to a situation which is out of the scope of this work. 



Let us finally point out that the special case /i = leads to a Langevin equation of the form 

dx ^ , , 



(B.14) 



which is again non-acceptable. In this case, however, detailed balance does not suffices to completely fix the properties 
of L{x,t). Although all the odd cumulants are zero, the even ones remain undetermined, allowing therefore for non- 
Gaussian noises. 



For odd variables, we have to replace e by —1 in eqs. ( |B.l| ) and (B.2). In this case, the first term in eq. ( p3.2| ) 



vanishes. Therefore, the homogeneous system of equations is already triangular, showing that at = for i > 3, unless 
the determinant is zero. This happens only if 

Vy^O (B.15) 
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I n th is case, nothing can be said about the set of cumulants since in view of eq. (B.15), in eqs^ (B.l) and 

( B.2 ) all the coefficients vanish. The corresponding Langevin equation also has the form given in eq. ( B.14 ), again 
non-acceptable. 
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APPENDIX: FIGURE CAPTIONS 



Fig.l: Alkemade's diode. Thick lines represent plates and wires of metal 1 while thin lines represent plates and wires 
of metal 2. The system enclosed by the frame corresponds to the diode and it is kept at constant temperature T. 
The figure on the right hand side represents a schematic view of the system. Note that the junction of metals 1 and 
2 give rise to a potential barrier. 

Fig. 2: Alkemade's diode in parallel with a condenser. To avoid the junction in the wires, the plates of the condenser 
are of different metals. We study the charge fluctuations q in the condenser of capacity C. 
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